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Abstract 

In this work we found the new class of exact stationary solutions for 
2D-Euler equations. Unlike of already known solutions, the new one con- 
tain complex singularities. We consider as complex, point singularities 
which have the vector field index greater than one. For example, the dipole 
singularity is complex because its index is equal to two. We present in ex- 
plicit form a large class of exact localized stationary solutions for 2D-Euler 
equations with the singularity which index is equal to three. The obtained 
solutions are expressed in terms of elementary functions. These solutions 
represent complex singularity point surrounded by vortex satellites struc- 
ture. We discuss also motion equation of singularities and conditions for 
singularity point stationarity which provides the stationarity of complex 
vortex configuration. 

1 Introduction 

The importance of exact solutions for 2D-Euler equations is well known. To- 
day, the list of exact solutions is quite impressive. Without pretending to be 
exhaustive we will mention only some of them. First of all, this are classical 



solutions with smooth vorticity, such as Rankine and Kirchhoff type vortices 
(see, for example, the standard references pQ- [3] ), elliptical Moore and Suff- 
man vortices. [2], Lamb dipole [T] and Stuart vortex pattern [3]. It is interest- 
ing, that these classical solutions are still topical (see for example recent works 
[SJ,[5]). The generalization of these solutions are the models of different coher- 
ent structures, vortex patches, and vortex crystals (see, for example, [7J-[I3] and 
references therein) , which are well observed in numerical and laboratory exper- 
iments (see, for example . [T4] - [2"T] ) . Let us note that Stuart solution is based 
on the equation of Liouville type for stream function. Others classes of exact 
solutions rest upon the equation of Sinh-Poisson type for the stream function 
[2"2"] ~ |24j . Some exact solutions in Lagrange coordinates are given on the work 
[25] . It should be noted the interesting class of Kida -Neu vortices [26] . [27] ■ 
Also there is a lot of solutions with singular distribution of vorticity (see. [2] 
), and solutions, which contain a smooth part of vorticity field as well as point 
singularities [35]- [32], which are usually forming symmetrical configurations. 
Numerous publications deal with particular cases of point vortices (main ref- 
erences can be found in [2], [3], [33] and in works [31]- [35J ). It is known, that 
point vortices generate integrable dynamical systems as well as dynamical sys- 
tems with chaotical behavior (see, for example, [3J, [33], [37], [35] )■ The case of 
a couple of point vortices, which vorticities are of same magnitude but opposite 
signs arouses particular interest. Such point vortices couples are usually called 
point dipoles. Many works are dedicated to the dynamics and statistic of point 
dipoles of this type (see, for example. [3] . [36] . [39]- [41].) Such point dipoles at- 
tract interest, especially in the theory of equations Charney-Hasegava-Mima for 
atmosphere and plasma. It means solutions of modon type [42] and its models 
obtained with the help of dipole structures (see, for example. |43j- [47]) . Such 
a point dipole has two elliptical singular points (FigfTJ). However, there are 
also 2D-point dipoles, which have different structure. They have one complex 
singular point which has the vector field index |J| = 2. This type of point 
dipoles appears naturally as a consequence of expansion of stream function 'J 
in multipoles moments, similar to usual electrodynamics. More precisely, we 
suppose, that the vorticity uj is localized in the restricted neighborhood G of 
the point r o, which belongs to 2D-plane. Let us denote the characteristic 
scale of domain G by a, i.e.||G|| < a. Then the stream function 



Here, the first term is obviously point vortex with the vorticity strength 
Fo : F = J ui(~ct)cfit . The second term is the dipole ^ d : 
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,we can expand, as usual, in multipoles moments : 
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Figure 1: Couple of point vortices with vorticity of opposite signs, which is 
considered as point dipole. 

with dipole moment = j w(a)acTa. Further, we will call this dipole 

by dipole singularity to avoid any confusion with point dipole which is composed 
of two point vortices. It is easy to see the difference between them when calcu- 
lating the vorticity. Apparently, the vorticity of vortices couple is proportional 
to the difference of two 6— functions, whereas the vorticity of dipole singularity 
is proportional to derivative of S— function: 

to ~ S (x) = lim — [S(x + h) — S(x)], 
However, if the parameter h is considered as minor, but final, 

lo ~ -\5(x + h)- 6(x)], 
h 

then such a couple of closely located point vortices with the opposite vorticity 
may be considered as a approximative model of dipole singularity. 

The important role of dipole singularities in dynamic of 2D-Euler equation 
was mentioned in the paper [48] . It is shown in this paper, that point dipole 
singularities are moving singularities of 2D- Euler equation which are form- 
ing together with point vortices the hamiltonian system of singulaties motion 
equation. This dynamical system has three independent integrals of motion in 
involution of Kirchhoff type. This fact means the complete integrability of the 
problem of motion of one dipole singularity and of one point vortex. Correspond- 
ing exact solutions for the plane case without boundaries are given in the paper 
|49] . As it was noted earlier, the dipole singularity is more complex than point 
vortex since the index of its vector field is equal to two ( unlike simple singular- 
ities, which index is equal to ±1). In this paper we will continue the studies of 
solutions for 2D-Euler equation with complex singularities, which were started 
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in [35] > Also, we discuss in details some questions which previously were 

presented briefly. Further, we demonstrate, that more complex multipole singu- 
larities, generally speaking, are not compatible with the dynamic of 2D- Euler 
equation. But in this work we show that 2D-Euler equation has a new class 
of exact stationary solutions with complex singular point, which index is equal 
to three. These solutions are found in explicit form and expressed in terms of 
elementary functions. Obtained solutions describe localized vortex structures, 
in which complex singular point is surrounded by vortices satellites. In addition 
we discuss the equation of motion for singularities and we give the sufficient 
conditions of immobility for singular points; without this equation one can not 
guarantee the stationarity of solution. 



2 Dipol singularities in 2D-Euler equation 

We use 2D Euler equation in form of Poisson brackets for the vorticity u> and 
the stream function * : 

u) = -A* (3) 
{A*,*} = 0, (4) 



dt 

or in the form: 

<9A* dA* dA* n 

here 1?- is fluid velocity. Poisson brackets has the usual form: 

{A B}-e dA dB - dA dB dA dB 
dxi dxk dx dy dy dx 

£ik~ is single anti symmetrical tensor: 

£12 = l,£2i = — l)£n = £22 = 0; 
^ = (x,y) = (xi,x 2 ). Velocity field f = (V x ,V y ) = (V U V 2 ) i 



is written m 



the form : 



5* 

Vi = £ lfc - — . (6) 
ox k 

As a stationary solutions for the Euler equation one use often the anzatz: 

A* = /(*), (7) 
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where /( V I / ) is arbitrary enough differential function of For example, for 
the Lamb solution pQ, f($?) is chosen as a linear function, for the Stuart solution 
[4 /(*) = exp(-*) and for solutions [22] /(*) = -sinh*. 

In solutions of type (JTJ) the vorticity lo is smooth enough function. There 
are also stationary solutions in which the vorticity has a smooth part as well as 
one dimensional immobile singularities of point vortex type ( see, for example, 



In order to find stationary solutions for Eulcr equation (UJ with singularities 
of point vortex type it was proposed in the work [32] the anzatz: 



For the simplest case exact stationary solutions were found independently 
and in a different way in papers [31] , |32j . In this work we study a more general 
anzatz. Let us suppose that the vorticity uj = — A 'I' can be presented in the 
form: 



may depend on time. Then it is supposed that vorticity is composed of 
smooth part and singular part as well, which is a generalized function. To 
begin, let us consider physical sense of development ([8]). First term in right part 
is obviously a smooth part of vorticity field. Second group of terms is vorticity, 
which is generated by the set of point vortices N, with stream function '3/ at: 



The third group of terms in development ([8]) with first derivatives of 8 - func- 
tion describes the vorticity which is generated by group of dipole singularities 



[28]-[3j] ). 



N 



a* = /(*) - r a <5(^ - ~^ a ) 



N M 





M: 
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Actually, the let us apply the operator 




(9) 



to Laplace equation: 



— Aln|^-^| 2 = <5(^-^) 
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It is obvious that: 



£a *'!■-(*£)'(*-*'>. ,10, 

That is why the vorticity which is generated by the third group of terms has 
the stream function $m : 



M 



ln|^-^| 



(11) 



i.e. that is the sum M stream functions of dipole singularities in the form: 



X — X f 



which are in the points xp and have dipole moments 
field of dipole singularity has evidently the form: 



(12) 



The velocity 



V D = — 



D 



dy 



2ir \~$ — Xq\ 



2y0^) 



V. 



D 



&V D _ 1 
dx 2tt|^-^ | 2 



D, 



2x0^) 



The presence of complementary sources of vorticity in form of derivatives 
of 8 - function, itself, is not forbidden in ((HJ , if they are compatible with Euler 
equation ^ . We have to remind that the singularity part of vorticity like every 
generalized function with point support is composed of S— function sum and 
its derivatives only. However, in the next chapter we show that the condition 
of compatibility of development (j8]) with Euler equation is not trivial and en- 
genders important restrictions for (j8]). Formulae which one needs to work with 
derivatives of generalized functions are well known, but for reader's convenience 
we gave them in the Appendix. Before studying the general expansion of vor- 
ticity ((8|) let us examine some particular cases. We begin with trivial solution, 
which is describing one point vortex. Let us remind, how from point of view of 
generalized functions theory, point vortex satisfies formally Euler equation. We 
will substitute vorticity and stream function of point vortex which are in point 
io = 0, in Poisson brackets. Then we obtain: 



2n(x 2 + y 2 ) 



[x8(x)}6 (y) 



r 2 



2n(x 2 + y 2 ) 



[yS(y)]5(x) (13) 



This Poisson brackets is equal to zero since the terms in brackets are equal 
to zero from the generalized function theory point of view. The fact that Pois- 
son brackets is equal to zero (fT3|) is physically interpreted as absence of self- 
interaction in point vortex. If the Poisson bracket does not vanish for the so- 
lution with one singularity on the plane, then, typically, a self acceleration of 



G 



this singularity appears . This effect is considered unacceptable from a physical 
point of view and the corresponding solutions must be rejected. 

Let us examine in the same way one dipole singularity, which is in the point 
ito and has the dipole moment D. 

The Poisson brackets ((H) for one dipole singularity gets the form: 

{ ^ )A ^} = ^--£S^^ - ? )-__U^^ - Xo ). (14) 

Let us show that Poisson brackets (fiH) is equal to zero. We substitute in 
(|14[) derivatives in explicit form and calculate it in details: 

J^^5(a^-a^o) = DiS xx (x - x )6(y - y ) + D 2 5 x (x - x )S' y (y - y ), (15) 



Q v, Q f t a 

^^D^^(^-^o) = D 1 S x (x - x )5 y (y - y ) + D 2 S{x - x a )S yy (y - y ). (16) 
&Z D _ 1 Di[(y - yo) 2 -(x- x ) 2 } - 2D 2 (x - x )(y - y ) ^ 



dx 2vr [t-~#o\ 
dVp = 1 D 2 [(x - x Q ) 2 -(y- y ) 2 } - 2D Y (x - x )(y - y ) 
dy 2tt |^~^o| 4 
As a result we obtain: 



(18) 



2n{^ D ,A^ D } 



(D 2 + D 2 )(y-y ) 



(D 2 + D 2 )(x-x Q ) 2 
|^-^o| 4 



j 5 v (y - yo) 



5 x {x - x Q ) 



S x (x-x )- (19) 
S y (y - yo)+ 



DJMty - yo) 2 ~(x~ x ) 2 } - 2D 2 {x - x ){y - m ) 
|^-^o| 4 



S(x - x )S (y ~ y )- 



DiD 2 [{x - x ) 2 -(y- y ) 2 } - 2D 2 (x - x )(y - y Q ) » 

— °(y - yo)°xx\ x ~ x o)- 

\x - ld \ 

Let us use formulae from the Appendix. Action of first derivative of 8— function 
on usual functions is given by formula (|95|) . It follows from this formula that 
the first and second terms in formula (fl9l) are zero, since they contain zeros of 
the following form: 



7 



(y - yo)$(y - 2/0); (y - yo) 2 \y- Vo =o S y (y - y Q ); 

(x - x )S(x - x ); (x - x ) 2 \ X - Xo =o S x (x - x ). 

Besides, in third and fourth terms the term with the factor D|, turns into 
zero because it contains zeros (x — xq)S(x — xq), and the term with the factor 
D\ turns into zero in the same way, because it contains zero (y — ya)S(y — yo)- 
In addition, it is evident that the terms: D\D2[{x — xq) 2 S(x — xo)]5 yy (y — yo) 
and 

-D1-D2K?/ — yo) 2 d{y — yo)]fi X x( x ~ x o) are equal to zero. Consequently the 
brackets (fH?)) is equal to: 



2tt{* d , A$d} = Dl ° 2{V . V f s{x - x )S' n (y - y Q )- (20) 

D x D 2 {x - Xq) 2 » 

— — _> , 4 °(y - yo)s xx (x - x ) 

\x — ¥q\ 

Now, for computing this equation we use the formula ([97)1 , which is needed 
to apply the second derivative of 6— function in the commutator (|20[) . It is 
obvious that all the terms turn into zero, excluding terms without derivatives 
of 5— function, which are mutually eliminated: 



Thereby we proved that there is no self-interaction in dipole singularity and, 
consequently, it is the exact stationary solution of Euler equation. We can 
also understand the absence of self-interaction in dipole singularity basing on 
simple physical considerations. Actually, from (FigJ5]) one can see, that due to 
symmetry of stream line configuration, the flux of impulse which is flowing in 
singularity through arbitrary section S is exactly equal to impulse flux which 
is flowing out from singularity through the same symmetrical section S. That 
means that the force which is acting on singularity is equal to zero and the 
singularity does not move. 

Naturally the question arises: do the higher multipole singularities, for ex- 
ample, quadruple, self-interaction? To answer this question wc have to calculate 
Poisson brackets {^f, A^} for quadruple singularity. The quadruple singularity 
has the stream function: 

= ~D ilia ^-4— In | a* - ^o| 2 , (21) 



and the vorticity: 



in 1 2 dxi 1 dxi 2 



u= k D "m^ s &-^ (22) 
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Figure 2: Dipole singularity with the index of vector field equal to 2. 

When computing the Poisson brackets {'J, A'f} using the formula (l98l) we 
obtain the following result: 



which shows that there is self-interaction in the general case in quadruple 
singularity; (exception to the rule is special choice of coefficients: Di 2 = 0, 
or Dn + D22 = 0). The self-interaction means that the given singularity is 
not physical. The similar result is obtained for multipoles of higher orders. 
Therefore, in general case the multipoles singularities of orders higher than 
dipole are not compatible with Euler equation. 

3 Singularities motion and stationarity condi- 
tions. 

According to the work [48] let us examine the question in which cases vortic- 
ity expansion (j8|) is compatible with Euler equation (HJ) and which restrictions 
arise from it for stationary and not stationary cases. We consider singularities 
coordinates ~af a ,!ctp and all coefficients r a ,Z)f for general case as depending 
on time. Further, we will distinguish two cases: /(W) = 0, ^ 0. Let us 
first examine the first case /(\&) = 0. In this case the smooth part of vorticity 
is absent and the stream function is determined by singularities only 



{*D a ,A*z, 2 } 



1 £>i 2 (£>ii + £>2 2 ) 
37r (a; 2 + y 2 ) 3 



8&), 



(23) 



(24) 



and derivatives of vorticity obviously have the form : 
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N 



N 



sa* _ _ ^ di- ^ _ ^ ^ ^-^°) d< 
at ~ ^ dt 1 j ^ a* dt 

9 2 (5(^ - ^) dzf 



a = l 



f d_p[ 8S& - ^) f n0 

^ dt dXi ^ 1 



/3=1 



0=1 



dxidxk dt ' 



(25) 



Vn 



dXr> 



N 



= -^E r 



a <9£ (^-^ a ) 
dXr> 



M 



dx ' i 



(26) 



Where V p - is fluid velocity © , which can be written down in the following 
form: 



^ T Q In [2 - ^"| 2 + jT(Df ^-)ln |^ - 



Ct = l 



/3=1 



(27) 



We substitute (f2"5"j) and (j2"6")l in Euler equation©, and we obtain: 



AT r 



£ 

a=l 



dt 



9a;,- 



dt 



(28) 



£ 

/3=1 



d£>f gg(g - j^) ^ d 2 £(^ - afe) dxf «^(]^-Ji^) 



dt eta; 



dxidxk dt 



dxi dxp 



Since equations for different singularities ~2? Q (t), x^(t) must go to zero in- 
dependently, we obtain: 



iv 

E 



di- ^ _ _^ Q) _ ra ^-^) d< + y ra ^-^») 

dt dt p dx p 



0, (29) 



^ d£>f tM(l> - ^) cP dH{^ - -P) dxj | v ^ dHtf - gf) =Q 
^-^ dt 9a;, ' dxidxk dt p 1 dxidx p 

(30) 

Every equation (|29|) . (|30p contains also singularities of different orders (8— functions 
and derivatives of S— function). To satisfy each of equation (|21)1) . (15!J|) the fac- 
tors of singularities of different orders shall also go to zero independently. Let 
us begin with the simplest equation (|2TJ)) . The derivative of S— function acts 
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on velocity field according to formula (|95l) . Because of non compressibil- 
ity condition -o- 2 - = 0, the third term in formula (|29|) does not give additional 



contributions to singularities with S— functions. That is why when factors of 
5— functions are going to zero this gives the well known equation: 

i.e. conservation of vorticity strength. Terms with the first derivatives from 
S— function give: 



N 

\ " 



-l^^^o, ( 31 , 



^ \dt p dx 

i.e. point vortices motion equations: 



I ' dx a \ 

("df " ypl « Q J =0,(a = l,2,...JV) (32) 

(In equations (|32|) there are no terms with self- interaction). Let us examine 
now equations for point dipoles (|30p . The third term in equation (|30p con- 
tains second derivatives from S— function. These derivatives act on velocity 
field V p (x,t) according to formula (|97|) in Appendix, and give the first ones 
and second ones as well. Since derivatives of different orders must be set equal 
to zero independently, the equation (|30[) splits into two groups of equations: 




Equations (|33)) . (PH) are satisfied, if all factors of different singularities turn 
into zero independently and the number of equations coincides with the number 
of variables. Terms of point dipole with self-interaction as this was shown in 
previous chapter are absent. As a result we obtain motion equation of point 
dipoles: 

dD? «dV p 



(35) 



dxl 
~df 



dt % dx 



V P \ M ,,(0 = 1,2,...M). (36) 



From (f3"3"| . (|3^|) follows that if all dipole moments are equal to zero, (Z?f = 
0, V/3), then equations (l33l) . (f34|) are absent and only the system of equation (|32l) 
remains for point vortices motion. And vice versa, one can see from equation 
(|3"Tj) that if all T a = 0, Va, then the system (|3"Tj) is absent and only motion 
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equations of point dipoles remain. For the general case the equation system 
p2|) , ([35l) . ([36l) . where the velocity has the form ([27]). describes motion of point 
vortices N as well as motion of point dipoles M. This equation system was 
obtained for the first time in [48) . where, in particular, was shown, that it can 
be written down in the Hamiltonian form. Further, to write down the equation 
system of singularities motion, we denote coordinates of point vortices as 
(where a— is the number of point vortex a = (1,2,. ..TV) , index i takes the 
values i — 1,2, i.e. denotes the coordinates of point vortex x\ = x,X2 = y, 
the subscript v means the vortex coordinates). In similar manner we denote 
the coordinates of dipole singularities as x (where ft— is the number of dipole 
singularity /3 = 1,2, ...M). Then, motion equation of singularities (f32|) . (|35l) . (|36|) . 
taking into account formula (|27|) for fluid velocity, takes the form of : 



dxi 



dt 



~£ik 



N 



^7 ( X vk x Zk) | 



7#" 



(37) 



M p 
7T 

0=1 



"Ik 



2{ x vl x dl)( X vk X dk) 



x v x d 



X ,, X J 



= 0, 



dx 



di 



dt 



-£ik 



N 



r Q {xP dk x" k ) 



a=l 



2tt 



X X , 



(38) 



M 

£ 

7^/3 



Vk 



2{x 



c dl)( x dk 



X j X ■ 



-$P _ 



d0f 



/V 



- nP 



ik 



X d 



~^ta 



2( x dfc ^SfcX^dm ^Sro) 



x d 



|8 

dm 
x d 



(39) 



M 7 
7#£ 



''dm/ 



2Sml (x^ k 



v dk) 



x d 



It: 



25mk(Xdl X dl) 



8(x 



L di 



){x 



' d.k 



x d 





dm 



dm ' 



X j X j 



x d x d 

Equations (|37f , (|38)l describe the motion of point vortices and dipole singu- 
larities as well and equation (f39f describes time evolution of dipole moment. As 
this is shown in work [48] these equations have the Hamiltonian form: 
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dr a BH 



dt - 9 Dl 
dD? OH 



~£ik- 



Here Hamiltonian H has the form: 



c=l J /3=l c*=l,0=l 



V- z -"ml x dm ■ l dm) lJ l \ x dl ^dl) ly m ly rn\ x d x d) 



/9-i7=l 



Equation system (|4"0)) has the conservation laws of the type of Kirchhoff's 
generalized integrals i.e. the conservation laws related to motion equation in- 
variance under translation and rotation of coordinates system: 



N M 



h = Y< r *<i - Yl D i = const - ( 42 ) 

a=l /3=1 
N M 

h = Y T * x v2 ~ Y2 D 2 = con st. (43) 



a=l f3=l 
N M 

J = ^r tt (^«) 2 -2^0^) = const. (44) 

Q = l /3 = 1 

As in case of point vortices there are three independent motion integrals in 
the involution: H , J and if + /| . This means that the motion of one vortex 
and one point dipole are integrable. The exact non stationary solutions for 
plane case without boundaries are given in [49] . Naturally, the question arises, 
is it possible or not to add for the non stationary case in expansion more 
higher derivatives of S - function, i.e.multipoles of higher order. Multipoles of 
higher order, than dipole produce two kinds of difficulties. First of all, as it was 
shown in sectoin 2, that, generally speaking, singularities of this kind have self- 
interaction. From dynamical point of view, the substitution of high multipoles 
in Euler equation in accordance with the formula (|96p Appendix A, engenders 
overdetermined equation system. Hence, for the non stationary case, vorticity 
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expansion © is compatible with Euler equation with /(\f r ) = 0, il the following 
sufficient conditions are satisfied: 

1. Multipole moments starting from the quadruple one are equal to zero. 

2. Singularities motion obeys the equations (|4"0)) . 

Let us now examine stationary case, when /(W) ^ 0, ^jffl ^ 0. In this 
case the stream function can depend on time since singularities coordinates and 
dipoles moments are function on time. The substitution of vorticity © in Euler 
equation ((SJ) gives for the smooth part of stream function the equation: 



dt h w dt Q ys dt J h \ dx * dt ^ dt J 

dm _ Q 

The first term is related to explicit dependence of stream function on time, 
the second one and the third one are related to singularities motion. The fourth 
term in (I45[) is related to dependence of dipole moments on time. We have to 
add to the equation (|45p equations for singularities parts of vorticity field which 
were already given (formulae (|40j) ) . In the case — 0, the equation (l4"5j) is 

absent and only singularities motion equations remain. The sufficient condition 
of stationarity, i.e. the (|45|) goes to zero when ^ consists in the following: 

1. Stream function does not depend explicitly on time. 

2. All singularities do not move. 

3. All dipole moments D± are not function of time. 

All terms which contain the velocity V p in equation (|28|) form obviously the 
Poisson brackets {\E', A^}. That is why the sufficient condition of immobility 
for all singularities and of stationarity for all dipole moments is the condition 
when Poisson brackets goes to zero: 

{*,A*} = 0, (46) 

This means that all factors of all independent singularities in Poisson brack- 
ets go to zero (HBl . 



A dDf 



dt 



4 Exact stationary solutions with complex sin- 
gularities 

Now we consider the problem of exact stationary solutions of 2D-Euler equation, 
when d ^jy ^ 0. Further, it is easy to consider the Poisson brackets 

9* dAV dA* _ 

dx dy dy dx 
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as dimensionless. Let us choose the anzatz 



in the form: 



A* = cxp + ^nT Q S(lt - lt ) - Anl^ -^6(lt - ~a? c 



) (48) 



when function f(*&) is chosen in the same way as in the Stuart work [4] ; we 
suppose that coefficients T and D are constant, and the coordinate Xq is not 
depending on time. (For the sake of simplicity we can choose Xq — 0). n— is 
positive integer number. By means of evident rescaling: 



— > 
D 



1' 



(49) 



equation 
ted): 



is reduced to the more simple equation (the primes were omit 

9 - (50) 



A* = cxp(-'J) + 4irn5( it) - 4TrD-^S(lt). 

(Let us note, that when D =0 the equation ([5U1) has the solutions found in 
[3T],[35]). First of all, let us find exact solutions for the equation (f5TH) . and then 
we prove, that they are exact stationary solutions of 2D-Euler equation (|4T|) . 
Now we can look for the solutions of equation ([50)) in the Liouville form: 



ln8- 



u (z) 



(i + Hz)\y 



(51) 



when u (z) = , is the unknown for the moment function of complex 



variable z = x + iy and u(z)— is primitive function. 
Direct calculation of A^ (|5~Tj) gives: 



A* = 



u'(z) 



A In 



u (z) 



(52) 



(l + \u(z)\y 

it is important to note that the formula (|52p is valid for arbitrary analytical 
function u (z) independently of its singularities structure. 

We substitute the formula (|52|) into equation (|50|) and obtain the equation 



for the function 



u (z) 
- A In 



i 2 > d 

u (z) = &%n5(Je) - 4nD-g=jf8(lt). 



(53) 



It is easy to see that the equation ([53)) is satisfied, if we choose the function 

2 



u (z) 



in the form: 



u (z) 



1 -A 9 1 I I 2 

— 7t- exp JJ — — r In \ z\ 

\z\ 2n M 



(54) 
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Indeed: 



In 



u (z) 



-Tl\sx\z\ + 



^JL\n\z\ 2 



d 



(55) 



The first term in (|55j) gives the Green function of Laplace equation: 

Aln|z| 2 = 4ir5(l?) (56) 

and describes the point vortex. The second term in (|55|) is a result of appli- 
cation of the operator ~t) to the equation (|5rJ|) and describes the point dipole. 
Let us introduce the complex dipole moment : 

D = D l +iD 2 . (57) 
Then the dipole operator ~D -^f can be written in the complex form: 



^ dz ^ dz ' 



(58) 



(When D 7 z - denote complex conjugated values). 
The function (|54[) can be written down in the form: 



u(z)\ 2 = ex p (d± + D^j (In z + In z) = ± exp 1 exp (§) . 

(59) 

From formula (|59[) follows, that functions u (z) can be chosen in the form: 

u'(z) = -Lexp(^). (60) 

In the point z = 0, the function u (z) (|60l) has essential singular point, which 
joins the pole of order n. Now let us find the primitive function u n {z) : 

u n (z) = j exp \jA ^. (61) 
Using the new variable w = ~ , we obtain: 

u„H = -^rr / W {n - 2) exp WW. (62) 



From the formula (|62j) we can see that the primitive function is elementary 
function only with n > 2. This particular case is examined in this work (others 
cases will be considered separately). 

Integration by parts in the formula (f62|) with n > 2, gives: 

1 (D 

Dz n-2 eX P ^ z 

where the polynomial P„_2(z) has the form: 



u n {z) = - n ^„_ 2 exp ( — ) P n - 2 (z); (63) 
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P n -2(Z) 



l-(n-2)(-) + (n-2)(n-3) - 



(64) 



+(-l)- 3 (n - 2)! (-i)"" 3 + (-l)- 2 (n - 2)! (J) 



n-2 



As a result, |u„(z) has the form 

1 



\u n {z)\ 



U|2(«-2) 



D 1 x + D 2 y 

ex P I 2 „2 ,„2 



(65) 



The primitive function %(z) (|63p . like the function u (z) (|60[) has in z = 
the essential singular point, which joins the pole of order (n — 2). In the real 



form 



u (z) 



, has obviously the following form 
2 1 



u (z) 



(x 2 + y 2 y 



D x x + D 2 y 

eX P 2 2 i 2 

a; + y 



(66) 



Consequently, the essential singular point describes in the complex form the 
singularities of point dipole kind, while the pole describes the point vortex, since 
the expression (1661) generates following terms in stream function (|5ip : 



In 



Hence, the exact solution of the equation ([5U|) is given by the formula (l5"Tj) , 
where u(z) is defined by the expression (|6"3")l , while u (z)— by the formula 
Now we can prove, that the obtained solution turns into zero the Poisson brack- 
ets (l47|) . At first, let us calculate the velocity field. Derivatives ^ and ^ have 



the form: 



9* _ 2 d 
dx l+\u{z)\ 2 dx 

d^ _ 2 d 
dy l + \u(z)\ 2 dy 



u(z)u(z), 



Using the formulae: 

d_ 

dx 



d_ + d_ d_ 

dz dz dy 



u{z)u(z) 



d d 



dz dz 



we obtain a more convenient formula for derivatives: 

2 



d^ 

dx 

a* 

dy 



1 



Kz)r 

2i 



_du du 
u- — h u— 

dz dz 



_du 

dz 



du 
dz 
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Taking into account the formula (|60[) . after simple algebraic transformations 
we obtain the expression for components of velocity field: 



dx 

d^ 

dy 



2 exp 



Dz+Dz 



(l + \u(z)\')\D\'\z\ 
2ie xp(^±^) 



^(DFP n - 2 +Dz*P n -i 



^(DrP n -2-Dz 4 P n - 3 



(67) 



(68) 



(Let us remind, that n > 2). Now we show that function (j5Tj) . (|63l) is an 
exact solution of Poisson brackets (|47|) . For that we substitute the expression 
for vorticity (f48|) and derivatives (|67l) . f[68]) into Poisson brackets (l47|) . First of 
all we examine the simplest case n = 2. In this case the polynomial P„_2 = 1 
and derivatives ([67| . (|68l) take the simple form: 



dx 

a* 

<9y 



Dz+Dz 



4 exp 
(1 + |u| ) \D\ \z\ 

Aexp (nj±Dz) 

77—, L„| 2 | , 4 (^ 2 - 2^y - £ 22/ 2 ), 
(1 + \u\ ) \D\ \z\ 

H z )l = 77772 ex P I 772 I • 



1^1 \ |*| , 

We write the Poisson brackets (|47|) in the explicit form: 

{vEf, A*} = inn 



<3* <9* / ' 

— (y) - — %)<5 (x) 



(69) 

(70) 
(71) 

(72) 



-47T 



1 dx 2 dy 



5 (x)S (y)- 



-47T 



3* » 3* 
£ 2 — S(x)5 (y)-D 1 —6(y)S (x) 



It is obvious, that all the terms in the first brackets (|T12|) are equal to zero 
because they contain this kind of zeros: 

x 2 S(x),xS(x),x 2 S (x); y 2 S(y),yS(y),y 2 6 (y). 

In the second brackets one part of terms is also equal to zero, but there are 
dangerous terms of this type: xyS (x)5 (y). However, these terms are part of 
second brackets in the following combination: 



'9* 9* N 

Dl ~8x~ ~ D2 ~dy~ I 6 ' {X}6 ' (y) = [■•'K- D i D * x V + D i D l xy)f (x)S' (y) = 0, 
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i.e. are reciprocally cancelled. (Here brackets [• • • ] denote the common 
factor). Now we consider the last brackets in (|72|) . In this brackets also, one 
part of terms turns into zero at once, but there are dangerous terms of this kind: 
y 2 S(x)5 (y) and x 2 8(y)5 (x). These dangerous terms are part of the brackets 
(jT2"j) in the following combination: 



■}[D 2 D iy 2 S(x)5 (y) - D 1 D 2 x 2 S(y)S (x)] 



(73) 



(Here brackets [■ • • ] denote the common factor). Now we use the formula 
(|97[) . From this formula one can see, that dangerous terms have the form: 



Dl D 2 5(x)6(y)-^y 2 - Dl D 2 6(x)5(y)-^x 2 



= 



(74) 



and are cancelled in commutator (|72[) . Others terms are obviously zero. 
Consequently, the Poisson brackets turn into zero for all singularities. Accord- 
ing to the results of chapter 3, this guarantees that singularities do not move 
and that the dipole moment i5 is conserved. It is proved that the obtained so- 
lution of equation ([50]) . is exact, stationary solution of 2D-Euler equation (j4"7]) 
with n = 2. let us consider now the general case n >- 2. In this case velocities 
([68]) contain the polynomials P n _ 2 (z) i P n 



(j67|) . ([68)1 contain the polynomials P„_ 2 (z) i P n - 2 {z) ([64]) . 

It is clear now that the additional powers z or z in these polynomials gener- 
ate in Poisson brackets zero terms only. The dangerous terms coincide only with 
the first term in the polynomial P n -2(z), i.e. unit. But these terms correspond 
to the case n — 2 , and are already been considered. Hence, we prove that 
formulae in (1511) . together with the function u n (z) (|63[) . give exact stationary 
solution of 2D-Euler equation with n > 2. In explicit form this solution has the 
form : 



ln8-lnu(z) + 2 ln(l + \u(z)\ 2 ) 



(75) 



ln8 + nln(a; 2 + y 2 ) -2 



D x x + D 2 y 
x 2 + y 2 



+ 



+2 In 



\Pn-i{z)X 



\Df{x 2 



y2)(n-2) 



exp 



D\x + D 2 y 

x 2 + y 2 



From equation ([75]) it follows that the essential singularity splits up into 
singularities of point vortex and point dipoles types. However, as we will see 
later, the fusion of these singularities leads to a singular point with more complex 
geometry with a vector field index equal to three. This can be interpreted as 
the sum of the indexes of the point vortex and the point dipole. 
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5 Examples of vortex structures with complex 
singularities 



1) First of all let us examine the simplest case, when n — 2. In this case the 
function 1*2(2) ((6"3")) has the form: 



u 2 {z) 



1 



D 6XP V : 



D 



The function |m2(z)| 2 is given by (|71l) . u 2 (^) has the form: 



u 2 (z) = [x'+y ) exp 2 



x 2 + y 2 



As a result we find the stream function ([FIT) 



ln8 + 21n(, 2 + , 2 )-2 J ^ + ^ - 



+2 In 



1 + ^2 Z n2 ex P ( 2 ^2 T ~2 



(76) 



(77) 



(78) 




Figure 3: The simplest vortex structure with the index of vector field equal to 
3 (n=2). Complex singularity is a result of the fusion of essential singular point 
(dipole singularity) with the pole (point vortex). One can see the presence of two 
hyperbolic points A and B, external and internal correspondingly. Separatrix of 
these points connect these hyperbolic points with the central complex singular 
point. Straight lines indicate exceptional directions of central singular point. 

Farther, we can consider that D\ = D^- The vortex structure which is de- 
scribed by the stream function ((75)1 . is presented on (Fig (3)), with D\ = D2 = 1. 
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It is clear, that with great values | af| 2 — > oo, the stream function (|75|) coincides 
asymptotically with the stream function of point vortex with negative vorticity. 
From (FigJ3J) one can see that inside of external closed stream line there is a vor- 
tex structure with the a non trivial topology of stream line. In the center there 
is a complex singular point, which results from the coincidence of singularities 
of point vortex and point dipole types. Furthermore, one can see the presence of 
two hyperbolic points, one outside A and another inside B . Separatrix of these 
points link the hyperbolic points with the central complex singular point. We 
examine more in details the structure of complex singular point. The general 
theory of such kind of singularities is stated in qualitative theory of ordinary 
differential equations (see, for example, [50], [SI] )■ According to this theory, 
first of all, we need to choose exceptional directions of the singular point. This 
are directions of tangents, following which the infinite number of integral curves 
go inside of the singular point and outside of it. One can see from (FigJ3]), that 
there are four such exceptional directions which are designated on (FigJ3]) by 
direct lines. Integral curves which are inside of exceptional lines form sectors. In 
our case, between the lines there are only four elliptical sectors (see, for exam- 
ple, [51] ). According to the general theory, the index J(0) of complex singular 
point is given by Bendixson's formula (see, for example, |50j): 

J{$) = \{2 + n e -n h ), (79) 

where n e — is the number of elliptic sectors and nu— is the number of hy- 
perbolic sectors. In our case nh — 0, n e — 4. That is why 

J(0) = 3. (80) 

The condition (|80p means, that the complex singular point is structurally 
stable, because the necessary and sufficient condition of the structural stability 
for complex singular point on plane is, that its index J(0) ^ 0. (see, for example, 
[51]). Note, that in case of dipole singularity exceptional directions coincide 
with vector direction of dipole moment while dipole index Jjj — 2, i.e. point 
dipole is the structurally stable singularity. Index of complex singular point can 
be found otherwise without using of the general theory. For this, let us cut out 
by circles singular points as it is shown on (Figf?]) . Then, in the obtained multi 
connected domain the index of vector field is equal to zero. I.e.: 

J* + + J s = 0, (81) 

i 

where J s — is the index of outside circle S, J2i Ji~ IS tne sum of indices of all 
internal simple singular points. This means, that the index of vortex structure, 
which is surrounded by contour S, is equal to: 

J A + Jb + J(0) = 1. (82) 

Since index of hyperbolic points A and B is equal to (—1), then the equation 
(l82|) . gives index of complex singular point J(0) = +3. 
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B 














A 





Figure 4: Contour, which is used to find the index of complex singular point. 



2) Now, let us examine the case n = 3 (dipole plus pole of order n = 3). In 
this case the polynomial P n - 2 {z) is not trivial: P\(z) = 1 — j-. 



Function 1 1x3(2)! has the form : 



or in the real form: 

\M z )\ 2 = 



\D\ 



exp 



Dz + Dz 



1 - 



D 



y 2 ){Dl 



■ exp 2 



Dix + D 2 y 



y 



1 



(x 2 +y 2 ) (D 1 x + D 2 y) 



- 2 



{Dt + DD {Dl + DD 
Correspondingly the stream function ^ has the form: 



* = -ln8 + 31n(.T 2 + y 2 ) - 2 



D lX + D 2 y 
x 2 + y 2 



(83) 



(84) 



(85) 



+21n[l + |u 3 (*)f], 

where ^3(2) | 2 is given by the formula (|84|) . Stream lines picture is presented 
on (Figj5j) with D\ = D 2 = 1. One can see that, unlike previous case, the 
elliptical point appears in the solution. Previous internal hyperbolic point splits 
into two. Structure of central singular point does not change, its index is still 
equal to 3. 

Note, that outside separatrix remains the same. It surrounds all internal 
vortex structure, including new elliptical point around which appears the vortex 
with negative vorticity. The separatrix of previous hyperbolic internal point 



22 



-1 -0.5 0.5 1 1.5 X 



Figure 5: Vortex structure with n=3 (dipole and pole of order n=3). One can 
see, that the elliptic point appears. Internal singular hyperbolic point splits 
into two. Structure of central singular point does not change. Straight lines 
indicate, as earlier, exceptional directions of central singular point. Structure of 
the external separatrix does not change. Structure of internal separatrix which 
are now surrounding new elliptic points becomes more complex. 

gets more complex because now it surrounds another additional elliptical point. 
All separatrix link hyperbolic points either with each other or with central 
singularity. 

3) Now we consider the case n = 4 (dipole plus pole of the order n = 4). In 
this case the polynomial P n -2(z) has the form: 




(86) 



In the real form: 



\P2(z)f 



= 1-4 




1 + 2 



{x 2 + y 2 ) 
{Dl + DD 



) 



+ 



(87) 




+ 4 



(x 2 + y 2 ) 
(Dl + D 2 



) 



1 + 



(x 2 +y 2 ) 
(Dl + Dl 



) 



) 




(88) 



As a result for stream function we obtain the expression: 
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Figure 6: Vortex structure with n=4, "butterfly" (dipole and pole of order n=4). 
One can see, that two elliptic singular points appear with two vortices satellites. 
Besides, one new group of separatrix appears, whiich connects vortices satellites 
with the central singular point. Topological structure of separatrix of external 
and internal hyperbolic points does not change. 

9 = - In 8 + 41n(x 2 + y 2 ) - 2 + ° 2V + (89) 
+21n[l + \ui(z)\ 2 }. 

The image of stream lines is presented on (FigJHJ) , with D\ = D2 = 1. 
(Butterfly). First of all one can see that two singular elliptical points appear 
and around them two satellites vortices. We can see also that there are four 
hyperbolic points. Separatrix structure becomes more complex. There are three 
groups of separatrix. First separatrix is of the same kind that outside separatrix 
in all previous cases. The second one is the same as the separatrix of internal 
hyperbolic point in first case. So a new group of separatrix appears which links 
vortices satellites (elliptical singular points) with the central singular point; the 
structure of the last one does not change from tropological point of view. 

4). Let us examine poles of higher order :n = 5, i n — 6. Correspondingly 
polynomials P n -2(z) have the form: 

P 3W = 1 _ 3 (J) +6 (£) 2 _ 6 (i.) 3 , m 

ftW _,_4(*) +U (*)"_ M (»)\ M (£)\ (91) 

Stream lines with n = 5 are presented on (FigJT]), with Di = D2 = 1. In 
this case three elliptical singular points appear (three vortices satellites) and five 
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Figure 7: Vortex structure with n=5. In this case three elliptic singular points 
appear (three vortices satellites) and 5 hyperbolic points. 

hyperbolic singular points. In case of n = 6, four elliptical singular points appear 
and six hyperbolic singular points. With the increasing of number n there are 
always n- hyperbolic singular points and n — 2 elliptical points. Central singular 
point conserves four exceptional directions, i.e. its index remains equal to 3. 
Appeared vortex structure has symmetry relative to square diagonal. Diagonal 
pass always by central singular points and by opposite singular point, which is 
hyperbolic with n = 2k, and elliptical with n = 2k + 1. The obtained vortex 
structure has the form of necklace composed of vortices satellites except the low 
sector, which always has hyperbolic singular point, linked by separatrix with 
central singularity. 

6 Discussions and conclusions 

In this work we want to call attention to the fact that there are exact solutions 
of 2D-Euler equation which contain point singularities more complex, than sin- 
gularities which are usually considered as typical ones for 2D-Euler equation. 
Such complex singularities can be non stationary |48) . or stationary as well. Let 
us remind, that complex singularities are defined as singular points of vector 
field, which have the index |J| > 2 . The simplest singularity of this type is the 
dipole one with the index J — 2. Point vortices and dipole singularities form 
a set of moving singularities in 2D-Euler equation, which dynamics is Hamil- 
tonian [48 . With /(VP) = 0, for general case, moving singularities can be only 
point vortices and point dipoles. It means, that its index can not exceed two. 
The reason for this, as it was shown earlier, that there are self-interaction of 
multipoles and overdetermination of their motion equations. 

Now let us examine the case when in expansion ([SJ the function f(*f?) ^ 0. 
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Consider more specifically the anzatz © in the form: 

A* = /(*) - T 5(^ - lt ) - ^^S(^ - ^o). (92) 

It is not difficult to see, that with /('J) = 0, there is no stationary solutions 
for 2D-Euler equation because {A\l/, \&} ^ 0. If the function f(^) ^ and is 
chosen in the Stuart's form: 



/(*) = exp(-*), (93) 

the situation change substantially, what can be seen from results of this work. 
Presence of smooth part of vorticity field in the equation (|92|) gives exact solu- 
tions with the more complex singularity of index 3. As it is shown in this work, 
the singularity of vector field of index 3 can be interpreted in complex form as 
a fusion in function u (z) (|60|) of the pole, which corresponds to point vortex, 
with essential singular point which corresponds to point dipole. In real form, as 
it can be seen from the equation (1751) , for the final stream function the essential 
singularity splits up into singularities of point vortex and point dipoles types. 
However, the fusion of these singularities leads to a singular point with more 
complex geometry with a vector field index equal to three. This can be inter- 
preted as the sum of the indexes of the point vortex and the point dipole. Exact 
localized solutions obtained in this work describe vortex structure of the com- 
plex form, where the singular point is surrounded be vortices satellites. With 
increasing of the number n the vortices satellites have tendency to form sym- 
metrical necklaces. The existence of exact solutions with complex singularities 
itself is an important fact that is why in this work we contented ourself with con- 
sideration of simplest class of exact solutions expressed in elementary functions. 
We did not deal with questions of the construction of more complex solutions 
expressed by special functions and with important questions of stability of vor- 
tex configurations with complex singular points. All these questions must be 
examined separately and some of them will be studied in next works. 

7 Appendix. Action of derivatives of 5— function 
on velocity field 

For convenience we present in this Appendix some formulae used in this work. 
It is well known (see, for example, [52] ), that generalized functions are linear 
functionals, which are acting in the space of basic functions {<p(x)}. By def- 
inition, the derivative 5 (x) acts on differentiable function a(x) according to 
formula: 



S (x)a(x)(p(x)dx = — I S(x)a (Q)(p(x)dx — / S(x)a(0)ip (x)dx = (94) 



S(x)a (0)(p(x)dx + / 5 (x)a(0)<p(x)dx;\/<p(x). 
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This gives a well known formula: 

a(x)S' (x) = ~a (x)S(x) + a(0)S' (x) (95) 

Derivative of order (k) acts on infinitely diffcrentiablc function a(x), accord- 
ing to formula (|96p . which is obtained in a similar manner as a result of 
integration by parts: 

k 

a(x)S^(x) = J2(-iy +k Cia {k - j \0)S^(x). (96) 
2=0 

We need to use now particular cases: k = 2: 

a(x)5" (x) = a" (0)S(x) - 2a (0)6' (x) + a{0)s" (x), (97) 

and k = 3 : 



a{x)8" (x) = -a'"(0)S(x)+3a"(0)d'(x)-3a'(0)S"(x)+a(0)s" (x). (98) 

From formula (|96|) one can see that the derivative of order (fc) from 5— 
function engenders also all derivatives of lower orders including terms without 
derivatives, i.e. simply S— functions. Similar formulae are obtained also for the 
case of multiple variables. In particular, for two variables x^ , Xi 2 the second 
derivative acts according the formula: 



G> 2 

dx% ^ dx^ 



dxi± dxi 2 
da(x 



\-£=o 



50)- 



dxi 



860) 



dxi 



a(0) 



da0) 
dxi 2 

d 2 s0 



dxi x 



(99) 



0xi 1 &Xi 2 
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